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Multi-plasmon absorption in graphene 
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We show that graphene possesses a strong nonlinear optical response in the form of multi-plasmon 
absorption, with exciting implications in classical and quantum nonlinear optics. Specifically, we 
predict that graphene nano-ribbons can be used as saturable absorbers with low saturation intensity 
in the far-infrared and terahertz spectrum. Moreover, we predict that two-plasmon absorption and 
extreme localization of plasmon fields in graphene nano-disks can lead to a plasmon blockade effect, 
in which a single quantized plasmon strongly suppresses the possibility of exciting a second plasmon. 
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The field of nonlinear optics ranges from fundamental 
questions concerning light-matter interactions to excit¬ 
ing technological applications 0. However, usually very 
large field intensities are required to observe nonlinear 
effects. One is thus always looking for systems that will 
exhibit nonlinear phenomena at lower powers, with the 
ultimate limit being strong interactions between just two 
quanta of light [2]. One possibility to increase nonlinear 
effects is to use the strong localization and enhancement 
of electromagnetic fields in the form of surface plasmon 
excitations |3(. In that regards, we note that graphene 
dj has been demonstrated to support extremely localized 
plasmons HHI2. While optical nonlinearities in graphene 
have been studied by several authors jTTfl [2f| , here we pre¬ 
dict a novel nonlinear effect in the form of multi-plasmon 
absorption. We also show how this effect leads to sat¬ 
urable absorption in graphene nano-ribbons at low input 
powers in the far-infrared and terahertz spectrum. More¬ 
over, we predict that the extreme localization of plasmon 
fields in graphene nano-disks leads to such a strong two- 
plasmon absorption that it becomes nearly impossible to 
excite a second plasmon quanta in the system. This plas¬ 
mon blockade effect would cause the nano-disk to behave 
essentially like a quantum two-level system, which is ob¬ 
servable in its resonance fluorescence spectrum. 

Graphene is a two-dimensional hexagonal lattice of 
carbon atoms Q- The low-energy band structure of 
graphene is described by Dirac cones with the electron 
dispersion E n k = nhvF |k|, where vf = 10 6 m/s, and 
n = ±1 stands for the conduction (valence) band [25|. In 
its intrinsic form graphene is a zero-gap semiconductor; 
however, it can also be easily doped with free carriers 
and as such it supports plasmon modes At low 

frequencies, one can get a rather accurate description of 
these modes by using a simple Drude conductivity, 
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where Ep is the Fermi energy of graphene and 7 is a 
phenomenological damping rate that takes into account 
various decay channels like impurity or phonon scatter¬ 


ing 0. The resulting plasmon dispersion is given by 
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and we assume the average dielectric permittivity E r ~ 
2.5, which roughly corresponds to the case of graphene 
on SiC>2 substrate and air on top. This simple Drude 
model breaks down at large frequencies when plasmons 
become strongly damped by electron-hole excitations, 
which can be described by the Random Phase Approxi¬ 
mation (RPA) @- 0 . However, at low energies the Pauli 
principle blocks this decay channel and the plasmon is a 
long-lived excitation [1 - 121. The resulting plasmon wave¬ 
length A p = 2Tr/q p is about 100 times smaller than the 
free space wavelength A = 2i tc/uj, leading to the extreme 
localization of electromagnetic fields [?]. 

An intuitive explanation of the strong nonlinearities 
associated with plasmons emerges by considering the typ¬ 
ical doping levels in graphene. For an electron density of 
n = 10 12 cm~ 2 , the distance between two electrons is 
r e = l/yjrvn = 5.6 nm, and so to observe some kind of 
nonlinear phenomena we need to compete with an intrin¬ 
sic electric field of the order 
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This is significantly smaller than the characteristic field 
amplitude associated with nonlinear effects in atoms 0 , 
given by the field between an electron and proton at a 
distance of a Bohr radius as = 0.5 A. The electric field 
in that case is E at = 5 • 10 11 V/m, about 4 orders of 
magnitude larger than the field E e in graphene! 

To see what happens to plasmons at such a field 
strength E e , let us imagine a general case of plasma oscil¬ 
lations at frequency uj and wavevector q, which is accom¬ 
panied by an electric field E{ r, t) = i E p e l<i ' r ~ xult +c.c. in 
the plane of graphene. When the plasmon field is small, 
it will induce a surface charge density p v = q<T ^ E p that 
in turn creates an electric field E p = 2 _* , thus driving 

self-sustained charge density oscillations. However, at 
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FIG. 1: (a) Plasmon dispersion relation. Gray area denotes 

the regime where a single plasmon can excite an e-h pair. 
Red area denotes the regime where this process is forbidden 
but plasmons can decay via two-plasmon absorption. Sim¬ 
ilarly, the yellow area denotes the regime of three-plasmon 
absorption, (b) Two-plasmon damping rate 7 
for E p = E e , indicating that plasmons cannot oscillate when 
the plasmon field is equal to the intrinsic electric field. In¬ 
set shows the band structure of graphene and a two-plasmon 
absorption process. 


the field strength E p = E e , the induced charge density 
\p p \ = en/2 will be comparable to the initial charge den¬ 
sity en, and this simple picture of plasmons breaks down. 
We will in fact see that at this field E e there is a strong 
plasmon damping via the multi-plasmon absorption. 

To understand how this comes about, let us first look 
at the linear (single-plasmon) absorption. Assuming that 
the graphene plane is perpendicular to the z axis, the 
plasmon field can be described by the scalar potential 
<p(r,z,t) = |(/? p e lq ' r_9 l 2 le _lajt + c.c., where r = (x,y). 
The electric field is then given by E = — V</?, with the 
amplitude \E p \ = q\<p p \, for both in-plane and out-of- 
plane components. The interaction of the plasmon held 
with the electrons can be described by the Hamiltonian 
H p = — \ep v e lcl ' rj -c.c., where Tj is the position 
operator acting on the j th electron. To calculate linear 
absorption we can write the dissipated power as P = 

where 

rf,,/ 1 ) O-TT 

-^L = T \(n\H p \0)\ 2 6(E n -E 0 -nio) (4) 

is the Fermi golden rule for the probability that the sys¬ 
tem will absorb one plasmon quantum of energy fru. Here 
| n) is the many-body excited state of momentum fi,q and 
energy hio, with respect to the ground state 10), and we 
assume that the system is at zero temperature [26]. To 
quantify absorption we can look at the dissipation rate, 
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Here W = \ f d 2 rp(r)ip(r) is the total electrostatic en¬ 
ergy stored in the system, which is given by W = 
\e r £o\E p \ 2 A/q , where A is the surface area of the 
graphene flake. Finally, the relevant figure of merit is 
the plasmon quality factor Q = w/ 7 . 


From Eq. 0 we see that the absorption process con¬ 
sists of a sum over all events where the plasmon can ex¬ 
cite a single e-h pair. Conservation of energy and mo¬ 
mentum require that Hui = -E n2k+q — E rl|k : however, the 
Pauli principle allows this process only above the thresh¬ 
old condition of hio > 2 Ep — hvpq (see the gray area 
in Figure Q] (a)). We can calculate the linear absorp¬ 
tion by usingthe basis of Dirac electron wavefunctions 
in graphene [261 ]. As an example, at energy hio = 1.5 Ef , 
the dissipation rate 7 ^/ 0 ; = 1.1 is so high (Q = 0.9) that 
the plasmon is not a well-defined excitation. Below this 
threshold the Pauli principle blocks the absorption pro¬ 
cess and the plasmon is a long-lived quasi-particle. How¬ 
ever, if we increase the plasmon held, higher-order (non¬ 
linear) absorption must be accounted for as well. 

We note that this simple calculation of the linear 
absorption gives the same result as from the RPA 
[2(|. Encouraged by this fact, we proceed to calcu¬ 
late the nonlinear, two-plasmon absorption by writing 
dw m 

P = 2 _, tj 2 Hu and using the Fermi golden rule for the 
probability that the system absorbs two plasmon quanta: 
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Alternatively one could perform third order expansion of 
the single particle density matrix including the screen¬ 
ing fields consitently in every order of the expansion [27 1 . 
Such a calculation yields additional terms that contribute 
to the dissipation, arising from higher-order screening 
correlations. However, at high fields, when the nonlin¬ 
ear absorption is large, the screening process will be less 
effective and the simple calculation <[ 6 |) should give a good 
estimate of the absorption. By evaluating expression © 
we get the two-plasmon absorption rate 
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where F ^ is a dimensionless function of plasmon fre¬ 
quency, which is given in the Supplemental Material [26 ] 
and shown in Figure |T] (b). It is straight forward to show 
that the Pauli principle allows the two-plasmon absorp¬ 
tion only above the threshold hio > Ep — hvpq (see the 
red area in Figure [T] (a)). Then if we look at the en¬ 
ergy hio = Ep , single plasmon can not excite an e-h 
pair but it can decay via two-plasmon absorption. At 
this energy F^ 2 ' 1 = 1.6 and we see that intrinsic field E e 
sets a natural scale for the nonlinear effects. Specifically 
at the held E p = E e the two-plasmon absorption rate 
7 ( 2 )/w = 1.6 is so large that plasmon is not a well de¬ 
fined excitation. In fact at this intrinsic held we expect 
that the whole perturbation theory should fall apart. In¬ 
deed at the same energy fno = Ef and held strength 
E p = E e we obtain the three-plasmon absorption rate 
yl 3 ) /u) = 0.7 (26|; which is about 40% of the two-plasmon 
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FIG. 2: (a) Transmission spectroscopy of a graphene nano¬ 

ribbon. (b) Absorption cross section normalized to the surface 
area of the ribbon. Solid line stands for the low pump I -C I g , 
and dashed line for the high pump intensity I s = 6 kW/cm 2 . 
The saturation of absorption is caused by two-plasmon ab¬ 
sorption inside of the ribbon. The ribbon width is D = 25 nm, 
doping n = 10 12 cm -2 , and resonance frequency ujq = 2-7TX28 
THz. 


absorption rate, signalling the breakdown of the pertur¬ 
bation theory. Moreover, at the same field strength but 
lower energy hoj = 0.6 Ep (in the yellow area in Figure 
Q] (a), where two-plasmon absorption is forbidden by the 
Pauli principle) we obtain the three-plasmon absorption 
rate /oj = 6.5 again showing that plasmon can not 
oscillate at the intrisic field E e . 

These effects could be easily observed in experiments 
involving extended graphene as the broadening of the 
plasmon linewidth with increasing field amplitude. How¬ 
ever, we will now show that much more exciting be¬ 
haviour can be seen by using graphene nano-structures to 
provide a resonance and field enhancement for the plas¬ 
mon modes. The response of graphene nano-structure 
can be described by the polarizability [H, El : 

a(i0) — D 3 -r;- j -t— pr. ( 8 ) 

a ' K_ 4n£ouoD ' * 

E r cr(u;) 

Particularly, in the case of a nano-ribbon, K = 16 and 
G = L/D , where D is the width and I > D is the 
length of the ribbon El- To describe the plasmon res¬ 
onance we can use the simple Drude model of the sur¬ 
face conductivity cr(uj) given by Eq. 0 - However, now 
we must include the total damping rate jt = 7 + 7 nl, 
which contains both the linear term ( 7 ) like impurity 
or phonon scattering, and the nonlinear term ( r y n i) like 
two- or three-plasmon absorption. Specifically, to pro¬ 
duce a resonance at frequency ujq, we require a ribbon of 
width: D = 477 A p « A p /2, where A p = 2n/q p is the plas¬ 
mon wavelength in the extended graphene given by Eq. 
0 - We are primarily interested in the regime Huiq = Ep 
where a single plasmon can not excite an e-h pair, but 
two-plasmon absorption is allowed. For a typical doping 
n = 10 12 cm -2 (Ep = 0.12 eV), the corresponding free- 
space wavelength is Ao = 27rc/wo = 10.6 pm, plasmon 
wavelength A p = 62 nm, and a ribbon width D = 25 nm. 

The absorption cross-section is given by a a (cu) = 
47 r^Ima(w) = JfA 2 (mo 2 ^+^ /4 , where we have in¬ 


troduced the radiative decay rate of the ribbon = 

16-7T 3 e r G D 3 , . 

3 —TC 


20. Ei To estimate impurity or phonon 


scattering we can use measurements of direct current mo- 
bilit y 7 = ^ 0- For typical graphene mobilities of 
p = 10 4 cm 2 /Vs 0, we have 7 = 8.6 THz, while 7 ,. = 6.1 
GHz, even for a ribbon of length L = 100 D. Since 
7 r "C 7 , graphene nano-structures primarily act as ab¬ 
sorbers, while the absorbed power can easily be detected 
via the reduction of light transmitted across the ribbon 
iEl (see Figure [2] (a)). 

Dissipated power is P = Icr a = 7 t W, where W is the 
total energy, and I = ^£qc\E\ 2 is the incident light in¬ 
tensity. We can then estimate the plasmon field inside of 
the ribbon by using the result for the extended graphene 
W = \e r £o\Ep\ 2 A/qp , where A = LD is the ribbon area. 

We obtain: |_E7 p | 2 = A ^ o _J^a + 7 a / 4 |E| 2 , where the total 
damping rate 7 * = 7 + 7 n i itself depends on the plasmon 
field through the nonlinear damping term 7 n i oc \E p \ 2 
given by Eq. 0 - Then by increasing the intensity, there 
is an increase in the total damping rate, and a decrease 
of the absorption cross section (see Figure[2](b)). In par¬ 
ticular, at intensity I s = 6 kW/cm 2 , we obtain 7 n i = 7 
on resonance, which reduces the absorption cross section 
by a factor of 2. This corresponds to an input power 
of only 2 mW for a laser focused to a diffraction-limited 
spot size (A/2) 2 , which would induce negligible heating 
of the graphene flake [26|. On the other hand, extended 
graphene can also be used as a saturable absorber, due to 
Pauli blocking, but the required intensities are 4 orders of 
magnitude higher, I s = 0.3 GW/cm 2 0. Moreover, by 
lowering the doping we can further reduce the saturation 
intensity of the ribbon, while the resonance frequency 
would fall into the terahertz spectrum. 

Especially interesting is the case of a nano-disk where 
we can localize the entire field to an extremely small vol¬ 
ume ss A p . The polarizability of a disk can also be de¬ 
scribed by the expression 0, if we now take D as the 
disk diameter, K = 12.5 and G = 0.65 0. To produce 
a resonance at energy Huiq = Ep and doping n = 10 12 
cm -2 , requires a disk diameter D = 20 nm. This yields a 
radiative decay rate of 7 r = 24 MHz, while the other pa¬ 
rameters are the same as in the case of a ribbon (above). 

We can now estimate the electric field amplitude as¬ 
sociated with a single quantized plasmon by writing 
W = fujjQ, and using the result for the extended graphene 
W = T;£ r £o\E p \ 2 A/q p: where A = 7tD 2 /4 is the disk area. 
This gives a remarkable result for the field amplitude, 



2 8 7T 4 hio 0 
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2 • 10 7 V/m. 
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In other words, the field of a single quantized plasmon is 
of the same order of magnitude as the intrinsic field E e \ 
Then, for two plasmons the damping rate due to two- 
plasmon absorption would be so high ( 7^/073 = 5.6) 
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FIG. 3: (a) Resonance fluorescence of a graphene nano-disk. 

(b) Quantized plasmon energy levels in a disk. Two-plasmon 
absorption induces a huge linewidth of the doubly-excited 
state 12 ), which makes it extremely difficult to populate this 
state and effectively turns the disk into a two-level system. 

(c) Absorption cross-section normalized to the disk surface. 
The solid line stands for the low pump field (fl <C 7 ), and the 
dashed line for the strong pump (fl = r )/\/2). (d) Scattered 
spectrum on resonance (to = too) in the strong pump regime 
(n = 47 ), clearly showing Rabi sidebands at too ± fi. (e) The 
second-order correlation function g^(t) of the scattered field 
for weak pump intensities, showing an anti-bunching effect 
(J 2 ' 1 ( 0 ) <C 1 characteristic of a two-level system. 


that the resonance peak would completely disappear, 
leading to a plasmon blockade effect. 

To quantify this effect we adopt a density matrix ap¬ 
proach, -Jj = j^[p,H\+C [p\. Here the system Hamil¬ 
tonian is given by H = fuo 0 N — 4 p(ae l “* + a + e~ lut ), 
where N = a + a is the plasmon number operator, and 

f 2 = 7 r is the Rabi frequency describing the in¬ 

teraction with the incident field of intensity I. The Liou- 
villian C = C W + C'' 2 ' > consists of two terms that charac¬ 
terize the linear and nonlinear dissipation, respectively, 
£(*> \p\ = ^(2 a N pa+ N - a+ N a N p - pa+ N a N ) 0, 
where 7W = 7 and 7/wo = 0 . 05 , so that 7U) <C 7^. Un¬ 
der conditions of weak external driving, the dynamics of 
the (infinite-dimensional) density matrix can effectively 
be truncated to a few-excitation manifold [ 3 l|. Specif¬ 
ically at fl < y( 2 ) the steady-state population of the 
excited state |2) (containing two plasmons) is extremely 
weakly populated, P22 = 1+ ( 27 ^) my , and disk effectively 
behaves as a two-level system. 

The absorption cross-section is then given by <r Q = 

7r 7/2 


^ImPro = If A 2 


)2 +7 2/ 4+ q 2/2 > which saturates at 
intensity I s = 50 kW/cm 2 (Q = 7/7/2). Like in the 
case of the ribbons one finds that the radiative damp¬ 
ing rate is negligible, 7 r -C 7, and the disk primarily 
acts as an absorber. However, the weakly scattered light 
will now show interesting spectral properties, since it is 
well-known that a two-level system behaves as a strong 
frequency mixer [ 311 ]. To substantiate this, let us look 


at the power spectrum of the scattered light, S(to s ) = 
ARe f^°(a + (0)a(t))e i ^ at dt. By using the quantum re¬ 
gression theorem [3l| we obtain the steady-state two- 

time correlation function on resonance, (a + (0)a(f)) = 
\ (e -7 */ 2 + i e -37t/4 e -*nt + 1 e -3 7 t/4 el nt^ ewhere 

we have assumed the strong-pump regime SI 7 / 4 , but 
also S7 -C 7 ^ so that the disk still behaves as a two- 
level system. The scattered spectrum consists of three 
peaks; one at the laser frequency to 0 and two at the Rabi 
sidebands too ± B. In Figure [3] (d) we plot the case of 
S 2 = 4y, while tt/'yC 2 ' 1 = 0.04 so that the disk still be¬ 
haves as a two-level system up to an excellent approxi¬ 
mation. A peculiar feature of this result is that the sys¬ 
tem produces frequency mixing (characteristic of disper¬ 
sive nonlinearities) starting only from the two-plasmon 
absorptive nonlinearity. 

The two-level nature of the system is especially re¬ 
flected in the second-order correlation function g^ (t) = 
(a + ( 0 )a + (t)a(t)a( 0 ))/(a + ( 0 )a( 0 )) 2 , which describes the 
probability of detecting a second scattered photon at 
time t given a detection event at to = 0. In steady-state, 
low pump regime Q <C 7 , it is straightforward to show 
that </ 2 )(0) = (l + 7^ 2 Vt) 2 = 7 • 10“ 5 . Thus, the disk 
exhibits an almost perfect anti-bunching effect character¬ 
istic of an ideal two-level system (to compare, g^ ( 0 ) = 1 
in the absence of nonlinearities). The temporal duration 
of this anti-bunching dip around zero time delay t = 0 is 
given by ~ I /7 (in particular, within a two-level approx¬ 
imation g( 2 \t) = 1 + e -7 * — 2 e d / 2 jjnj), after which 
the system returns to a stationary value of g^ 2 \t) = 1 as 
illustrated in Figure [3] (e). The sub-Poissonian nature of 
the scattered light, </ 2 )( 0 ) < g^ 2 \t), is a distinctly non- 
classical feature that reflects the inability of a two-level 
system to emit two excitations simultaneously [32| . 

While we have focused primarily on the absorptive 
nonlinearity, graphene will show also dispersive nonlin- 
earities such as the Kerr effect [l|. Such a mechanism 
would in principle result in a shift in the resonance 
peak, in addition to the saturation of absorption in the 
graphene nano-ribbon. On the other hand, we do not 
expect that this would change our results in the case of 
a nano-disk, since this dispersive nonlinearity also leads 
to a plasmon blockade effect [2l|. Interestingly, the es¬ 
timated value of the dispersive interaction strength [ 21 ] 
is two orders of magnitude lower than the two-plasmon 
absorption rate. 

In conclusion, we have shown that graphene possesses 
a strong nonlinear response in the form of multi-plasmon 
absorption, which leads to the saturation of absorption in 
graphene nano-ribbons, and the plasmon blockade effect 
in graphene nano-disks. 
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Supplemental Material for ”Multi-plasmon absorption in graphene” 


The electron-plasmon interaction can be described by a Hamiltonian: H p = —^etp p (u>,q)p+e lult +c.c., where ip p is 
the plasmon potential at frequency w and wavevector q , while p q = e _ * qr3 is the Fourier transform of the many- 

particle density operator p( r) = ]>^(5(r — r^) (Sl| . To calculate iV-plasmon absorption we can write the dissipated 

j 


(JV) 


power as: P ^ ^ NTuo , and use the Fermi golden rule to find the probability of absorbing N plasmon quanta 


w n0 


dt = 1 -\(n\H p \0)\ 2 5(E n -E o -hio) 


(SI) 


dw, 


( 2 ) 

nO 


dt 


2ir 


E 


(n\H p \m)(m\H p \0) 


Em En huj 


S(E n — Eq — 2 Huj) 


(S2) 


dw. 


( 3 ) 

nO 


dt 


2tt 

T 


(n\H p \m) (m\H p \l) (l\H p \0) 


(Em — E 0 — 2hio)(Ei — Eq — hw) 

m.l 


S(E n — Eq — 3 tuo) 


(S3) 


Here | n) is the many-body excited state of momentum Nhq and energy NHlu , with respect to the ground state |0), 
and we assume that the system is at zero temperature. In the most simple description of plasmons, we factorize 
our many-body state into the product of single-particle states which evolve in the mean-field, self-consistent plasmon 
potential p p (w,q) = - 5 IHIHQ- 

Low energy single-particle states in graphene are described by Dirac wavefunctions VVik(r) = (r|nk) = 

—= ^ ^ e ,k r with linear dispersion E nk = nfiVF\k\. where k is the electron wavevector, n = ±1, and we need to 

include the two-valley and two-spin degeneracy of the ground state [S5| . Finally, the damping rate can be written as 
7 (JV) _ pW/w where W = \e r E 0 q\ip p (w,q)\ 2 A is the plasmon electrostatic energy, and it is straightforward to show 
that: 


vW .... E„ 2N ~ 2 


1 -=fW( w ) 


E e 


(S4) 


where \E p \ = q\<p p \ is the plasmon field, E e is the intrinsic electric field given in the main text, and F are 
dimensionless functions of plasmon frequency given by: 


4-7T 


F^\w) = - - - V |(n 2 k + q|e ?q ' r |nik)| 2 ( 5(£ , ri2k+q -F; ril k-^)/r ll k(l-/r l2 k+q) 

E r £oq A ^ 


niri2k 


(S5) 


F^(w) = 


( 00 V E F 71-2 V- 

k F A n^k 


E 


( 713 k + 2 q|e* q ' r |n 2 k + q)(n 2 k + q|e lqr |nik) 


^n2k+q E n ik JEj 


X 8(E n 3 k+2q E n 1 k 2fiw)/ rai k(l /n 3 k+ 2 q) 


(S6) 


F ( 3 ) (w) = 


( w ^ 

10 Ep 3n 2 ^ 

( 714 k + 3q e lq r 7 i 3 k + 2q)(n 3 k + 2q e lq r 7 i 2 k + q)( 7 i 2 k + q e lq ' r 7 iik) 

\v F qJ 

k 2 F 2 5 A ^ 

r niri4k 

n2n3 (-^n 3k+2q -^nik 2/kxj) (£' n2 i c +q E'nik fthl} 


X ^(F'n4k+3q E n ik 3/lw) i / ni k ( f /n4k+3q) 


(S7) 
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Here A is the surface area of the graphene flake, kp is the Fermi wavevector, Ep is the Fermi energy, and /„k = 
Q(Ep — .EVik) is the Fermi-Dirac distribution at zero temperature. Also we have explicitly written the expression 
(IS5l) to be evident that = 21m e(uj,q), where e(w, q) is the dielectric function calculated in the Random Phase 

Approximation fsii , [sij , while we have simplified expressions (IS6I) and (IS7I) by using the plasmon dispersion relation 

_ 2lT£ 0 £ r h 2 , 2 

9 — e*E F W ’ 

Note that even at the room temperature T = 300 K and doping n = 10 12 cm -2 , Ef/ kT ss 4.5 and our calculations 
based on the zero-temperature approximation will still be qualitatively valid. To be more specific, at energy Huiq = Ep 
where we were looking at two-plasmon absorption, we will also see a small contribution from the single-plasmon 
absorption. However, this effect will be suppressed by roughly e _A / fcT ps e -2 3 where A « Ep/2 is the energy gap 
required to enter the regime of single-plasmon absorption (see Figure 1 (a)). 

From a technical stand point, one might also worry about severe heating of the graphene flakes at high laser powers. 
For example, in the case of a graphene nano-disk, the saturation intensity is I s = 50 kW/cm 2 , with the corresponding 
absorption cross section on resonance: er a (u;o) ~ 0.25A (see Figure 3 (c)), which would induce heating of the disk at 
I s a a /A = 13 kW/cm 2 . However, note that graphene was reported to heat up at a rate R = 3.3 K cm 2 /kW when 
subjected to a high direct current flow, while sitting on a room temperature SiC >2 substrate which acted as a heat sink 
jsij ]. Thus, even in the continuous wave regime a laser would increase the temperature of the disk by only AT ps 40 
K. Moreover in the case of ribbons, the saturation intensity, and the induced heating is even smaller. 
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